4. Image Enhancement in the 
Frequency Domain 


* Chapter 4 (Image Enhancement in the 
Frequency Domain) 
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FIGURE 4-1 The function at the bottom is the sum of the four functions above it. 
Fourier's idea in 1807 that periodic functions could be represented as a weighted sum 


of sines and cosines was met with skepticism. 
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Fourier Transform of a real function is generally complex 
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Some Properties of Fourier 
Transform 


Convolution: 
The convolution in 


Кажа = | fla) g(x a)da “тона болап the 
ж g(3) o Вада 2222 
f(x) g(x) F(u) Glu) 
Correlation: 
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2-D Fourier Transform 
F(f(x y) =F(u v) = (| f(x yep- /2л(их+ lady 4 Flu, | e"? 
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Discrete Fourier Transform (DFT) 
1-D DFT pair for a signal sequence x(n) of length N: 


1 м: 1 - jaxnk 
X9 = M xne N , k=0,....N-1 
AV 


IN-1 J2xnk 
xD =% Х(Юе У, n-O,...N- 1. 


k=O 
2-D DFT pair fora "n seguen x(n,,n,) of area NxN: 


Pk nk 
Xs, k) Ea Y xn, , kk -0...,,N- 1, 
Ir9n- 
Jaink nk 
Xn,n)- 5, У ха, kje ~ , nmc-0..N-1 
| 1 
ometimes dem in inverse, ot; appears in each ‹ 


X(O, 0) = A: = > xn, m) = х(п, DC component 


д =0 m, =0 


Example: 1-02 DFT 


X(W = Sane Ў” SUN, k=0,..., № 1 


If Қап) = { 2, 3, 4, 4}: 


F (O) ==> fime = [2+ 3+ 4+ 4] =3.25 


n=O 


Ба) = fine j2z n/4 


n=O 
=-12+ Be ^"? + 4e 7 + Де 772] =- 0.5+ /0.25 


F(2) =- 0.25 
F(3) =- 0.5- 30.25 


Example: 1-D IDFT 
xu) => Хада У, n=0,..., N- 1. 


ЕСО) >s F(E =З.25+ (- 0.5+ 10.25) + (- 0.25) + (- 0.5- 70.25) 
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—3.254 C 0.54. /О.25)е””? + (- 025) eF + (- 0.5- /0.25) e/* 7? 
—3.254 (- 0.54 j0.25)(0+ р + (- 0.25)(- 1+ jO + C 0.5- /О.25)(О- p 
=3 


f(2) = Е Юе?" 24/4 


=3.25+ (- О.5+ jO.25) e* + (- 0.25) е” + (- O.5- 10.25) е” 
=3.25+ (- O.5+ jO.25)(- 1+ JO) + C 0.25) + JO) + (- 0.5- jO.25)(- 1+ 70) 
= 


f(3) Ex Раве” Зк 


—3.254 (- 0.54. 10.25) е”?З"/? + (- 0.25) e/** + (- 0.5- 10.25) e2 
—=3.25+ (- 0.5+ /О.25)(О- р + (- 0.25)(- 1+ 70) + (- 0.5- jO25(0- D) в 
-4 


Example: 1-D IDFT 


N-1 КУРСА 
хі) =» Х(Юе ХУ, n—O,..., N- 1. 
k=O 


f(n) =(2, 3, 4 А) 

= F(O)(1+ JO, 1+ JO, 1+ JO, 1+ 70) 

+ F )(1+ JO, O+ Л, - 1+ JO, O- JD 
+ F(2)(1+ JO, - 1+ 70, 1+ JO, - 1+ 70) 
+ P(3)a+ JO, O- JL - 1+ JO O+ Л) 


1-0 16-POINT DFT BASIS VECTORS 
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Another 2-0 DFT Example 


run 


a 
b 


FIGURE 4.4 

(a) SEM image of 
a damaged 
integrated circuit. 
(b) Fourier 
spectrum of (a). 
(Original image 
courtesy of Dr. J. 
M. Hudak, 
Brockhouse 
Institute for 
Materials 
Research, 
McMaster 
University. 
Hamilton, 
Ontario, Canada.) 
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Fourier Spectra 


dynamic range of Fourier Spectra is usually much higher t 
cal display devices to display. For display, we often use: 
Du, v) -C*logt*|F(u V), where Cis a scaling constant. 


Useful Properties of 2-0 DFT 


Periodic extension in "space" and "frequency". 


Linearity: 

аҳ(п,љ) + bx(n,n)- aX(k,k)+ bX(, А) 
Circular Convolution: 
Multiplication: 

х(ц,Һ,)% x(n,nm) Х(К,Е)* X;Gs, I) 
Shift (Translation): - 22 Xmmk* mk) 
xn- m,n- 13) = X(K,k)@e Ы 
Periodic іп Frequency Domain 

ЖК, К) =Х(К+ NK) =Х + М) =X(k+ МЕ + М) 
Conjugate Symmetry: 
X(k,k)- X C k,- k) = X (М k, № Е) 


2-D Convolution and Correlation 
2-D Convolution 
н 1 M-1N-1 | | 
f(x у* H x y) “MN A. Ёт) Ңх- my- п) 
f(x уж Ах, у © Flu v H(u, v) 
f(x Ух, у e Ни Vx Hu v) 
2-D Correlation 
f(x Y Ix =. У f nxt m yen 
f(x Jo Kx y е F(u vyu v 
f (x y)g(x y e Flu v- Gu v 
f(x y)» f(x y) e F (u Fu 9) =|Р(ц |^, powerspectrum 


£ l т ами к ع آلا اعلا‎ 


(convolution) 
اوم‎ 


One of the two 


Correct original images 


A+C = 1 


Missing 


В-Ы-1------ м 


Properly extended (padded) image 


Result of filtering in the frequency domain without 
properly padding the input images 


Р Correct 


Result of filtering in the frequency domain with 
properly padded input images. 


a b 

c 
FIGURE 4.38 
Illustration of the 
need for function 
padding. 
(a) Result of 
performing 2-D 
convolution 
without padding. 
(b) Proper 
function padding. 
(c) Correct 
convolution 
result. 


f(x,y): AxB 
h(x,y): CxD 


pad zeros to both 


to make them PxQ 


= А + С – 1 


B+ р – 1 
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FIGURE 4.41 

(a) Image. 

(b) Template. 

(c) and 

(d) Padded 
images. 

(e) Correlation 
function displayed 
as an image. 

(Г) Horizontal 
profile line 
through the 
highest value in 
(e). showing the 
point at which the 
best match took 
place. 


Highest correlation 


Gray-level 
profile line 


Separability Property of 2-0 ОЕТ 


The 2-D DFT can be done by performing row- 
wise 1-D DFT(or FFT) followed by column-wise 
1-0 DFT (ог FFT), and vice versa. 
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Fast Fourier Transform (FFT) 
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Complexity of 
Fast Fourier Transform (FFT) 


“Бог an N-point transform, ОНТ requires 
about N(N-1) complex additions, and № 
complex multiplications. 

FFT requires about N log N complex 
additions, and 
(N/2) log N complex multiplications. 


‘For N=8, DFT requires 64 complex 
multiplications, while FFT only requires 12 
complex multiplications. 
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Enhancement via Frequency 
Domain Filtering 


xn. m) = X У хк. RMN kn k 


The filtering can be done in discrete frequency 
domain: 


УСА, AD = ХС, RO HC A AD 


° With different lowpass cutoffs or highpass 
cutoffs, different 

filtering effects can be achieved. 
° Ideal LPF in ѕрайај, фута is not achievable. 
Unless we 


have infinite length , otherwise suffer 22 
from Gibhe 


Ideal Lowpass Filtering in 
Frequency Domain 


ыз) H(u, v) 


FIGURE 4.10 (a) Perspective plot of an ideal lowpass filter transfer function. (b) Filter displayed - 
image. (c) Filter radial cross section. 3 
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FIGURE 4.11 (a) An image of size 500 x 500 pixels and (b) its Fourier spectrum. The 


superimposed circles have radii values of 5, 15, 30, 80, and 230, which enclose 92.0, 
04.6, 96.4, 98.0, and 99.5% of the image power, respectively. 


а 
1111111 


ааа = a d d с 


Ideal 
Low-pass 
Filtering 


=m 


ШШШ Tin 


ааааааа | aa a а а-а-а ач 
ab FIGURE 4.12 (a) Original image. (Ь)-(Ғ) Results of ideal lowpass filtering with cutoff 
са frequencies set at radii values of 5, 15, 30, 80, and 230, as shown in Fig. 4.11(b). The 
е Е power removed by these filters was 8, 5.4, 3.6, 2, and 0.5% of the total, respectively. 
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Butterworth Lownpass Filters 
1 


1+[ Du, 9/ D" 


Hu, v) = 


D(u, v) 


yD(u,v) is the distance from the (0,0) frequency 


FIGURE 4.14 (а) Perspective plot of a Butterworth lowpass filter transfer function. (b) Filter displayed as an 
image. (c) Filter radial cross sections of orders 1 through 4. | 


tial Domain Butterworth Low-Pass Filte 


abcd 


FIGURE 4.16 (a)-(d) Spatial representation of BLPFs of order 1, 2, 5, and 20, and corresponding gray-level 


profiles through the center of the filters (all filters have a cutoff frequency of 5). Note that ringing increases 
asa function of filter order. 
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b FIGURE 4.15 (a) Original image. (Ь)-(Ғ) Results of filtering with BLPFs of order 2, 
d with cutoff frequencies at radii of 5, 15, 30, 80, and 230, as shown in Fig. 4.1106) 
f Compare with Fig. 4.12. 
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Gaussian Lowpass Filters 
H(u V) —е D (u,v)/2o* =e D (и,№/21№ 
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FIGURE 4.17 (a) Perspective plot of а СІРЕ transfer function. (b) Filter displayed as an image. (ў | 
radial cross sections for various values of Dp. 
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FIGURE 4.19 

(a) Sample text of 
poor resolution 
(note broken 
characters in 
magnified view). 
(b) Result of 
filtering with a 
GLPF (broken 
character 
segments Were 
joined). 
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Historically, certain computer 
programs were written using 
only two digits rather than 
four to define the applicable 
year, Accordingly, the 


company's software may 
recognize a date using "00" 
as 1900 rather than the year 


Historically, certain computer 
programs were written usin 
only two digits rather than 


company's gens may 1 


High-Pass Filters H,(uV =1- H,(u v) 
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FIGURE 4.22 Top row: Perspective plot, image representation, and cross section of a typical ideal highpass 
filter. Middle and bottom rows: The same sequence for typical Butterworth and Gaussian highpass filters. 
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oatial Domain High-Pass Filters 
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- HGURE 4.24 Results of ideal highpass filtering the image in Fig. 4.11(a) with Dy = 15, 30, and 80, 
- respectively. Problems with ringing are quite evident in (a) and (b). 


1tterworth High-Pass Filtering 
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FIGURE 4.25 Results of highpass filtering the image in Fig. 4.11(a) using a BHPF of order 2 with Dy = 1h 
30, and 80, respectively. These results are much smoother than those obtained with an ILPF. | 


ussian High-Pass Filtering 
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FIGURE 4.26 Results of highpass filtering the image of Fig. 4.11(a) using а GHPF of order 2 with D, = 15 
30, and 80, respectively. Compare with Figs. 4.24 and 4.25. 
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FIGURE 4.30 

(a) A chest X-ray 
image. (b) Result 
of Butterworth 
highpass filtering. 
(c) Result of high- 
frequency 
emphasis filtering. 
(d) Result of 
performing 
histogram 
equalization on 
(c). (Original 
image courtesy 
Dr. Thomas 

R. Gest, Division 
of Anatomical 
Sciences, 
University of 
Michigan Medical 
School.) 
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Band-Reject Filters 
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FIGURE 5.15 From left to right, perspective plots of ideal, Butterworth (of order 1), and Gaussian bandit 
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Band-Reject Filters 


FIGURE 5.16 

(a) Image 
corrupted by 
sinusoidal noise. 
(b) Spectrum of (a). 
(c) Butterworth 
bandreject filter 
(white represents 
1). (d) Result of 
filtering. (Original 
image courtesy of 
NASA.) 


Homomorphic Filtering 


An image f(n, ) is created by "reflecting" the light 
from an object t at has been "illuminated" by some 
light source: 


The iudi dad add 74, L5 P па «БУ, 


varying component (cont Въ ће overall dynamic 
range), 


n,n) isa fast varying component (affecting the 
local contrast). 


We like to decrease the low-frequenc: 
and amplify the high frequencies | 
(simultaneous dynamic range 
compression and contrast enhance 
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Homomorphic Filtering 


Separating these two components and de- 
emphasize/emphasize each component: 


Inx(n,1,) —lni(n,n)-lnr(n, n) 


Take Fourier Transform 


Ки v) = (u, v) + Ки») 
Filter with a filter function H(u,v) that de- 


emphasize the low-frequency and emphasize 
the high-frequency component. 
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Homomorphic Filtering 
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Original & Homomorphic 
Processed Results 
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